*Column 1		
use "$dataDir/dataset_fires_harvest_daily_weather", clear

		local depvar "nshf"
		local cluslvl "id_mills"

foreach var in `depvar'{
	reghdfe `var' afterXtreatment temp_celsius rain rain2, abs(id_mills i.year#i.month) vce(cluster `cluslvl')

	estimates store `var'_treat_did_wc
		qui sum `var' if e(sample)==1 & year<2015
		local mean1=r(mean)
		estadd scalar mean1 `mean1' : `var'_treat_did_wc
		estadd local moyearfe "Yes"
		estadd local millfe "Yes"	
		estadd local cluslvl "Mill"
		estadd local wc "Yes"
		estadd local poisson "No"

	}	

*Column 2				
use "$dataDir/dataset_fires_nonharvest_daily_weather", clear

		local depvar "shf"
		local cluslvl "id_mills"

foreach var in `depvar'{

	reghdfe `var' afterXtreatment temp_celsius rain rain2, abs(id_mills i.year#i.month) vce(cluster `cluslvl')

	estimates store `var'_treat_did_wc_nh
		qui sum `var' if e(sample)==1 & year<2015
		local mean1=r(mean)
		estadd scalar mean1 `mean1' : `var'_treat_did_wc_nh
		estadd local moyearfe "Yes"
		estadd local millfe "Yes"	
		estadd local cluslvl "Mill"
		estadd local wc "Yes"
		estadd local poisson "No"
		
	}				

	*TABLE DAILY (placebo):
		esttab nshf_treat_did_wc shf_treat_did_wc_nh  using "$mainDir/tables/TabA2_falsification_test.tex", replace b(%12.5fc) se(%12.5fc) keep( afterXtreatment) varlabels(afterXtreatment  "After 2015 $\times$ Regulated")  label starl(* 0.1 ** 0.05 *** 0.01) stats(mean1 N moyearfe millfe wc cluslvl poisson, fmt(%9.3fc %12.0fc) label("Pre 2015 mean" "Obs." "Month by Year FE" "Mill FE" "Weather controls" "Cluster level" "Poisson")) nonotes